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An exact solution of the model of fully packed loops of two colors on a square lattice 
has recently been proposed by Dei Cont and Nienhuis using the coordinate Bethe Ansatz 
approach. We point out here a simpler alternative, in which the transfer matrix is directly 
identified as a product of R-matrices; this allows to apply the (nested) algebraic Bethe 
Ansatz, which leads to the same Bethe equations. We comment on some of the applications 
of this result. 
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1. Introduction 



Models of fluctuating loops play a key role in two-dimensional statistical physics, and 
a range of well-known models (Ising, Potts, percolation, 0(n), to name but a few) can 
be conveniently studied through their reformulations as loop models. Exact results about 
loop models have been produced by a variety of techniques, including the Coulomb gas, 
conformal field theory, the Bethe Ansatz, and stochastic Loewner evolution. 

In this note we study the fully packed two-color loop model on the square lattice 
(henceforth referred to as the FPL'^ model) from the point of view of the algebraic Bethe 
Ansatz. The FPL^ model was introduced in |jl|] as a generalization of the four-coloring 
model of the square lattice edges [|] . It is defined by assigning one of two colors (black or 
white) to each lattice edge, subject to the constraint that every vertex be incident to two 
black and two white edges. In this way, the black and white edges form fully packed loops 
which are given fugacities and depending on their color. 

The FPL^ model has attracted much interest over the last decade. Successive advances 
in the Coulomb gas technique have permitted to compute the central charge and the critical 
exponents for a number of special cases: the four-coloring model (nfe,n^) = (2,2) 
the dimer-loop model (nb^n^) = (2, 1) 0, and the equal- fugacity case ni, = [^. This 
eventually led to the solution for general values of ni, and ^j. An interesting special 
case is that of Hamiltonian walks, with {rib, = (0, 1) [Q. A generalization of the FPL^ 
model, obtained by giving the loops a bending rigidity, was solved in 0. It contains as a 
special case the so-called Flory model of protein melting ||^. 

All these Coulomb gas results are obtained by making certain reasonable, but non- 
rigorous, assumptions about the long wavelength behavior of an associated interface model. 
The resulting critical exponents are however believed to be exact, and they have been 
successfully tested against numerical Monte Carlo and transfer matrix results. 

To give the results obtained by the Coulomb gas a rigorous status, and to go beyond 
it and derive results which are not obtainable from a continuum approach, it is natural 
to turn to the methods of integrability. Following earlier work on the four-coloring model 
||10|| , Dei Cont and Nienhuis [jll| have very recently succeeded in finding a coordinate Bethe 
Ansatz for the equal-fugacity FPL^ model. In particular they have computed the exact 
partition function. Moreover, they have shown that when Ub 7^ n^, the FPL^ model is not 
integrable, in agreement with earlier expectations [Q. 
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However the Coordinate Bethe Ansatz is a rather complex technique, which, in order 
to be made fuUy rigorous, would require investigation of a large number of specific con- 
figurations. The goal of the present note is to present a simpler alternative, in which the 
transfer matrix of the equal- fugacity FPL^ model (henceforth we note n = ni, = riw) is 
identified with a product of trigonometric R-matrices of Uq{sl{4)) with n = —q — (and 
an appropriate twist). Applying the (nested) algebraic Bethe Ansatz allows us to recover 
the Bethe equations of Ref. in a straigthforward fashion. As a bonus we obtain the 
central charge and the critical exponents, which are found to agree with the non-rigorous 
results of Ref. . We also comment on a n ^ —n symmetry of some of the sectors of the 
transfer matrix. 

The paper is organized as follows. The FPL^ model is defined in Sec. 2. In Sec. 3 
we define its transfer matrix in terms of an R-matrix that adds four vertices and a twist 
matrix that takes care of the boundary conditions. We then show how these matrices are 
related to those of the affine quantum group Uq{sl{4)) with alternating fundamental and 
conjugate representations. The corresponding Bethe Ansatz equations are discussed in 
Sec. 4, and we reproduce in particular the eigenvalues of the transfer matrix Finally, 
in Sec. 5, we give the expressions of the central charge and conformal weights for n < 2, 
and we compare our results to those obtained for fully packed loops on the hexagonal 
lattice. 



2. The FPL^ model 

Following 1^, we reformulate the FPL^ model as a 24- vertex model. For each vertex 
of the square lattice, the four incident edges are decorated with arrows of two possible 
orientations (outgoing or ingoing with respect to the vertex) and two possible colors (black 
or white), subject to the constraint that each of the four possibilities be represented exactly 
once around every vertex. Clearly, following the arrows of a given color traces out an 
oriented loop of that color. Note that reversing the arrows along any one loop, and leaving 
all other arrows unchanged, leads to another allowed configuration. 



i 



CO 



00 



Fig. 1: The six types of vertices in the 24- vertex model (the remaining 
vertices are obtained by 7r/2 rotations) with their corresponding weights. 
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We then assign a weight w = WhW^ to each vertex which is the product of weights Wb 
and Wyj coming from the oriented black and white loops respectively. The black weight 



Wb = oj (resp. Wb 



when the black loop makes a right turn (resp. a left turn) at the 



concerned vertex, and Wb = 1 when the black loop goes straight. For the white weight we 
choose the opposite convention, that is with w and exchanged. The weights of the six 
types of vertices which are unrelated by 7r/2 rotations are given in Fig. 1. 

The FPL^ model with fugacity n for both colors of loops is recovered by summing 
independently over the orientations of all loops (black and white). An anticlockwise (resp. 
clockwise) black loop contributes (resp. uj~'^) to the fugacity, as it must turn four times 
more (resp. less) to the left than to the right. Thus, u is fixed by 



n 



(2.1) 



In order to apply the Bethe Ansatz it is important to specify the boundary conditions. 
In the following we shall specialize to the case where the square lattice is wrapped on a 
cylinder, i.e., with periodic boundary conditions across a horizontal row of 2L vertices 
Note that the argument leading to ( pT|) only works for contractible loops, i.e., 
loops that do not wrap around the periodic direction. To obtain the correct weighing also 
for non-contractible loops one introduces a vertical seam separating the first and the last 
vertex in each row |Tl[] . Horizontal edges cutting the seam are assigned an extra weight of 
a (resp. a~^) when covered by a left-pointing (resp. right-pointing) arrow; the convention 
does not depend on the color of the arrow. Clearly, non-contractible loops can only wind 
once, so a is fixed by 

(2.2) 



n 



a + a 



<?- 



Fig. 2: Parity convention for vertices and edges of the square lattice. Even 
(resp. odd) edges are shown in dashed (resp. solid) linestyle. Even (resp. odd) 
vertices are shown as dashed (resp. solid) circles. 
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With these boundary conditions, the FPL^ model contains three conserved quantities 



TlJI . To explain these, we shall adopt a convention for the parity of both the vertices and 
the edges, as shown in Fig. 2. The three components of the conserved charge which are 
conserved by the evolution along the cylinder are 

Q=[l\-[ , (2.3) 

\lJ \N^i + NobJ 

where A^... is the number of vertical edges of a given parity (e = even, o = odd) in the 
concerned row, b (black) or w (white) refers to the color of the arrow, and t (up) or | 
(down) to its orientation. The constant term has been added for convenience. Strictly 
speaking, it would make better sense to talk about conserved charges with respect to a 
parity convention for the columns that does not alternate from row to row. In this respect, 
the charges ( p.3| ) only commute with the transfer matrix that adds two rows at a time. 



3. The transfer matrix 



I 
\ 
I 
I 



Fig. 3: Convention for the labeling of arrow states on an even edge. 

Before going on we shall adopt a convention for labeling the arrow state i = 1, 2, 3, 4 
of each edge in the FPL^ model. This is shown in Fig. 3 for the case of an even edge; the 
convention for an odd edge is similar, but with all arrows reversed (i.e., i — 5 — z). 



Fig. 4: The FPL^ model R-matrix adds two even and two odd vertices as 
shown. The arrow indicates the transfer direction. The parity of vertices and 
edges follow the conventions of Fig. 2. The four dangling edges below (resp. 
above) the thin dashed line specify the in-state (resp. the out-state). 

We have seen above that from the point of view of the conserved charges, it is most 
natural to build up the lattice by adding two rows at a time. We have also assumed 
that the horizontal strip is of even width 2L. In order to define the row-to-row transfer 
matrix T, we first define a 256 x 256 matrix, which we denote, in analogy with integrable 
models, by R; it adds four vertices as shown in Fig. 4. Using the weights of Fig. 1, it is 
straightforward to write R explicitly, in the basis obtained by using the labeling of Fig. 3 
for the external lines and by taking the tensor product of the corresponding vector spaces. 

The transfer matrix that propagates the system in the upwards direction then reads 

T = tTaRaL---Ra2Ral{^~^<»^), (3.1) 

where the subscript a denotes the "auxiliary space" (the pair of horizontal lines) of di- 
mension 16 and the subscripts 1, 2, . . . , L correspond to the L pairs of vertical lines which 
form the "physical space". The twist (8) f2 is a matrix in the auxiliary space which 
takes care of the effect of the seam; explicitly, Q = diag(l/a, l/a,a,a) acts on the upper 
horizontal line whereas fl~^ acts on the lower line. 

We now introduce another R-matrix which is related to the affine quantum group 
Uq{sl{4:)), where q — —u)~^ (so that n = —9 — 9"^), and which we call R. It is schematically 
described by Fig. 5, in which the two representations □ and | of Uq{sl{A)) appear in an 
alternating fashion. Conventions are such that at a vertex where two lines intersect, the 
first factor in the tensor product refers to the leftmost line of the in-state, when seen along 
the transfer direction. 
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□ 



□ 

Fig. 5: R-matrix for two rows and two columns, which alternatingly carry 

the fundamental representation of [/g(s/(4)) and its conjugate. The arrow 
indicates the transfer direction. 

The four R-matrices which appear on Fig. 5 can be expressed as: 



X 



+ {qx ^ — q ^x)P\ 



□(g)n5 



Ro®|(a;) = {q^x - q-^x-^)P°^^ + {q^x-^ - q-^x)Pl^^ 
P-^®n{x) = {q^x - g-2a;-^)P^n + (g^x"^ - q-^x)P^^^ 
Rg®|(2;) = {qx - Q"^x"^)P^g + {qx~^ - g"^x)P^^g. 



(3.2a) 
(3.26) 

(3.2c) 

(3.2d) 



We have as usual denoted R = IIR, where 11 is the operator that permutes the two factors 
of the tensor product. The P are intertwining operators which can be computed using 
representation theory; the parameter x is the ratio of spectral parameters of the horizontal 
and vertical line. We shall give the explicit matrix representations of the R-matrices below, 
after fixing the values of x. 

In the present context, one can first define separately transfer matrices for even 
and odd rows (however, only their product will be directly related to the previous 
transfer matrix T). Indeed, we define Rg(x) = Rg(g,n(x/xn)Rg(g,g(x/X|) and Rn(x) = 
Ro®n(^/^n)r^8)|(^/^|)5 where the products are meant as in Fig. 5, and the spectral pa- 
rameters of vertical lines supposed to be fixed. Define next 



Tg(x) = trgR^i(x) ■ ■ ■R^2(x)Rgi(x)0-i 
Tu{x) = tVoRaLix) ■ ■ ■Ra2ix)Raiix)n 



(3.3a) 
(3.36) 



where, as before, the indices determine the spaces on which the matrices act. Note that Q 
and can be considered as the same element of the Cartan subalgebra of Uq{sl{4)), but 
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in fundamental and conjugate representations respectively. Due to the Yang-Baxter equa- 
tion, the T|(a;) and Ta{x) form an infinite family of commuting matrices; their product, 
the two-row transfer matrix T{x,y) = Ta{x)Tn{y) = Tn{y)Ta{x), is itself of the form 



T{x,y) = tTaHaLix^y) ■ ■■'Ra2{x,y)Rai{x,y){0. ^ (8) n) 



(3.4) 



where R(a:,y) = Ra(a;)Rg(|y). 

The claim is that the two R-matrices for two rows and two columns R and R that we 
have introduced are related. At this point we choose all horizontal lines (whether odd or 
even) to have the same spectral parameter, and similarly for all vertical lines, so that the 
ratio is constant and is: x = q~^. Note that at this special value, the matrices Rq^d and 
Rg(g,| become proportional to projectors onto the antisymmetric sub-representations (this 
enforces the fact that two loops of the same color cannot cross each other). 

Explicitly: (c = — q) 



( 



RD®q(g) = c 
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(3.5a) 



(3.56) 
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(3.5c) 



(3.5d) 



Note that we use the standard bases for □ and | (which are dual bases of each other). 
One can then check that 

Yi = c^URU-^ (3.6) 

where the constant takes cares of the extra factors in Eqs. ( |3.5| ), and U is a, diagonal 
matrix that fully factorizes as a tensor product over the four incoming lines: U = Uhn ® 
Uhn® Uyn ® with as a possible choice 



a, 
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(3.7) 



Consequently, the corresponding transfer matrices T and T are also similar up to a 
constant: 

T = c'^^UyTU-^ (3.8) 



where Uy is the tensor product of Uyn and Uya for all vertical lines. 
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Obviously, space does not permit us to reproduce the resulting 256 x 256 matrices R 
and R. Rather, Fig. 6 gives three examples of matrix elements. Note that our convention 
for the R-matrix (see Fig. 4) is to keep fixed indices for the horizontal and vertical lines. 
Thus, if the arrow configuration (coded as in Fig. 1) is P4P3P2P1 for the out-state (read 
from left to right when looking along the transfer direction), it is P2P1P4P3 for the in-state. 
The three examples in Fig. 6 then read explicitly: -R8i,i8 = + uj~^; -Rio3,91 = 
and i?239,i88 = uj"^ + The corresponding entries of R are found from ( p.6|) . 



4. Algebraic Bethe Ansatz 

The set of commuting transfer matrices Ta{x), T|(x) can be diagonalized using the 
so-called nested Bethe Ansatz. We shall not describe this procedure here and refer to 
l^,n,15] for details. The eigenstates are built by action of operators which depend on 



parameters that we call u\^\ with i = 1, 2, 3 and /c = 1, . . . , m^^\ on a reference eigenstate 



(highest weight state) which has only white arrows pointing up. These parameters satisfy 

(0 

equations, which, in our parameterization, are algebraic in the e^^^fc , where 7 is such that 
q = —e~'^^. Explicitly, call uj^'^^ the diagonal elements of the twist O in the fundamental 
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representation (here, u^'^^ = u^"^^ = 1/a, uj^^^ = uj^'^^ = a), with 11^=1 ~ ^^'^ define 



the functions Q^'\u) = [1^=1 sm-f{u - 4'^), i = 1, 2, 3. = = 1. Then the Bethe 
Ansatz equations read 



for l<^<3, 1<A;< m^^^^ The functions Z*-*-* depend on the representations and 
spectral parameters of the physical space (and on the twist); here, one easily computes 
f^^\u) = w(i)(sin7wsin7(w- 1))^, f^'\u) = wW(sin7(M + 1) sin7(w - 1))^ for z = 2,3, 
f^'^^u) = uj^'^\sm'y(u + 1) sin7w)^. 

The corresponding eigenvalues of T|(w) and Tn(u), in the parameterization x = 
_e»7("+i)^ are 

tu{u) = X'^^^u) + X^^\u) + X^^\u) + X^^\u) (4.2a) 
t^{u) =X(i)(w) +X(2)(t,) +1(3)(^) +xW(w) (4.26) 

whereX«(.) = ^^^^^^f^^\u) andX«(-) = g^;^^|S^/(-^V). 

We now choose w = 0, so that the X^^\ X^, X^^\ vanish (this, once agam, can 
be interpreted as a consequence of the requirement that two loops of the same color do 
not cross each other). Finally, we consider the two row- matrix T = TnTg. Its eigenvalue 
is obtained by taking the product of the remaining terms in Eqs. ( ^4.2| ); rewriting the X^'^^ 
as functions of Q^*-*, getting rid of the extra factors sin 7 which compensate the factors of 
c in Eq. (p.8|) , we find the eigenvalues of T to be 



g(2)(0)2 u(^) \Qi^){0)J Q(i)(0) Q(3)(l) ^ a;(2) V g(2)(o) Q(i)(-1) Q(3)(0) 



_,Q^ / Q(^)(-1)Q(3)(o) Q(2)(-1) / Q(i)(0) Q(3)(i) \ 

Q(2)(o)y Q(i)(o) Q(3)(l)+'' Q(2)(0) ^ Q(^H-l) Q(^KO) J ^ ^ 

where in the last line we have used the explicit expression of the twist. This is precisely the 

square of the expression found in for the one-row transfer matrix. The correspondence 



of notations is as follows: 

Uk = 2i{u'j}^ + 1/2) Vk = 2i{u'j^^ - 1/2) Wk = 2iu^p (4.4) 
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Furthermore, the Cartan subalgebra produces three conserved quantities; in funda- 
mental and conjugate representations, they have the foUowing expression: (basis of the 
dual of the root lattice) 



QlD 
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(4.5) 



Combining this with the correspondence given by Fig. 3, it is easy to identify them with 
the components of the charge Q of Eq. ( p.3| ). Their value is determined by the numbers 
m*^*^ of Bethe roots of kind i: each root of the kind i decreases by one the i^^ component 



of the charge, starting from the reference state which has Q 
Eq. (|2.3|), we deduce that 




Comparing with 



(1) 



m 



m 



(3) 



= N^^ + Nob 



(4.6) 



Note that if rrS^^ > but m^^^ = 0, only even arrows are modified compared to the 
reference state (i.e. all odd arrows are white pointing up); and similarly for m*-^-* > 0, 



m 



(1) 



and odd arrows. 



5. Results and conclusions 

We have found that the nested algebraic Bethe Ansatz can be used to solve the FPL^ 
model at Uh = Uw The latter is therefore identified with a standard integrable vertex 
model associated to Uq{sl{4)), and its transfer matrix embedded into an infinite set of 
commuting transfer matrices; and many results follow immediately. 

In particular, the long distance behavior of this type of models has been studied by 



many methods (see for instance [Il6| , p!7| , [18| ) . For |n| > 2 the spectrum has a gap and the 
correlation length is finite. In the following we focus instead on the critical regime |n| < 2, 
and we parameterize n = 2 cos 7. 

As our model is isotropic, we expect the largest eigenvalue of the transfer matrix to 
have the following asymptotic behavior |[T9|] 

TIC 

logtoiL) = -Lfo + — + ■ ■ ■ for 00, (5.1) 
bL 
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where c is the central charge of the infra-red conformal field theory. Here note that only 
fundamental representations (□ and g) are used, so we are dealing with a non-fused model. 
Assuming the usual form for the ground state, standard computations (see e.g. [^) lead, 
in this type of models, to the following form of c 

{w\C-^\w) (5.2) 



7r(7r — 7) 

where r and C are respectively the rank and the Cartan matrix of the underlying algebra, 
and is a vector with components Wg = \ log(c<;*'*''/c<;*'*"'"^'*) that parameterizes the twist. If 
we now specialize to and to our choice of boundary conditions: wi = = 0, W2 = —27, 
we obtain 

c = 3-12— ^ -, (5.3) 

7r(7r — 7) 

which coincides with what was found in . Note that this is not the central charge of the 
^^(^3) conformal field theory — the latter can also be obtained within the framework of 
this model, but with a different twist. 

One can also investigate the nature of excitations above the ground state. They are, of 
course, gapless and describe solitons associated to the three fundamental representations of 
A3 interacting with the standard S- matrices [ P^^pHj] , plus possible bound states in certain 
regimes of 7. In the infra-red limit the dispersion relation can be linearized and the 
corresponding low-lying excited states are related to the conformal weights of the 
aforementioned CFT via 

Tr(r — 24A 1 

\ogUL) = -Lfo+ ^ +••• for 00. (5.4) 

bL 

One can check that the weights thus obtained fit with the formulae of the Coulomb gas 
picture: 

1 1 

A„ = - (e| i^"^ |e - 2eo) + - (m| |m) (5.5) 

where K = |(1 — 7/7r)C, e (resp. m) is the electric (resp. magnetic) charge which belongs 
to the lattices of weights (resp. roots) of A3, and cq is the background charge, related to 
our twist by eo = -^w. This constitutes a confirmation of the results of 0. Incidentally, 
the fact that for rib 7^ "^w the quadratic form K appearing in the conformal weights as 
given in is generically not related to a Cartan matrix can be considered an indication 
of the non-integrability of the model. 
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We have made some numerical checks of the structure of the ground state and excited 
state described above. As it turned out, this picture was confirmed for ti > 0; however, for 
n < the ground state of the usual form (with, using the notations of Eq. (^^), real w^, 
Vk, Wk) is not the state corresponding to the dominant eigenvalue of the transfer matrix. 
Indeed, we have found a n ^ —n symmetry of the eigenvalue spectrum corresponding to 
the sectors where both L — m^'^^ and m^^'> -\-m^'^^ are even, cf. Eq. (|4.6D; this applies thus in 



particular to the ground state sector which has m^^^ = m^'^^ = m^^^ = L. This symmetry 
can be described in the Bethe Ansatz equations as the transformation of the Bethe roots 
(with the notations of Eq. ([4.4|)): 



7Wfc ^ -(tt - 7)ufc IVk -{ti - 'l)vk -fWk ^ -{tt - ^)wk + ITT, (5.6) 

where we recall that n = 2 cos 7, so that —n = 2cos(7r — 7). One can check that this 
transformation leaves the eigenvalues ( [4.3| ) invariant. In particular, we conclude that the 
"real" ground state eigenvalue to only depends on |n|, and identifies with the one described 
above only for n > 0. For n < 0, the "fake" ground state eigenvalue corresponds to a 
state very high above the real ground state (even the bulk part being different as L — > 00), 
whose only special property is that it is the analytic continuation of the n > ground state 
eigenvalue. At the moment, we do not have a satisfactory explanation of this phenomenon. 
The non-unitarity of the n < theory, or our boundary conditions (we can only consider 
the theory on a cylinder, but not on a torus due to the issue of winding loops) might play 
some role in it. 

[We remark parenthetically that the symmetry of the ground state sector holds true 
more generally for the Ub 7^ Uw FPL^ model, under the transformation (n^, n^) 
{—Ub, —n^). This is based on the observation that, with suitable periodic boundary condi- 
tions, all terms in the high-fugacity expansion of the partition function have Nb + N^ even, 
where A^^ (resp. N^) is the number of black (resp. white) loops. To see this, represent the 
dominant state at nb,nw —> 00 as an "ideal state" in the four-coloring picture with 
black (resp. white) loops being an alternation of colors 1 and 2 (resp. 3 and 4). Note that 
the high-fugacity expansion of the 1-2 (black) and 3-4 (white) loops (disregarding their 
orientation) can be obtained by only permuting the colors around the two other types of 
small loops (say, of types 1-3 and 2-4), and that these loops stay of length four. Examining 
all possible 1-2 and 3-4 loop environments of a plaquette occupied by the 1-3 and 2-4 
loops proves our statement.] 
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It is interesting to compare our results with the work of Reshetikhin [21| on fuUy- 
packed loops on the hexagonal lattice. After contracting the vertices of the hexagonal 
lattice two by two, so as to form a square lattice, this author identified the R-matrix with 
that of the integrable model Uq{sl{3)). The choice of spectral parameter, as in our case, 
makes R degenerate into a projection operator. However, there are important differences. 
First, in the underlying algebra is of course different and all horizontal and vertical lines 
carry the same representation □ of Uq{sl{3)), in contrast to the alternation of □ and g used 
here. Second, in there is no twist O in the auxiliary space, whence contractible and non- 
contractible loops carry respective weights of n and 2. In particular, for n < 2, the central 
charge is constant, c = 2. When n = 2, the continuum limit of the hexagonal-lattice loop 
model becomes a SU{3)k=i free field Wess-Zumino-Witten theory; this is a consequence 
of the Uq{sl{3)) identification and of the fact that only the fundamental representation is 
used. Likewise, the Uq{sl{4)) identification of the FPL^ model reported in the present work 
implies that the n = 2 case is a SU{4)k=i WZW theory in the continuum limit. Indeed, 
the four-coloring model was originally constructed by Read so as to have a SU{4)k=i 
symmetry. 
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